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PROGRAM 



Double Precision Signed Multiply 



TAPES 



ASCII Sources 090-000016 



ABSTRACT 

This routine multiplies two double precision, fixed 
point, two's complement numbers to form a quadruple 

precision t two's complement product. 
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1 .1 Memory 

IK or larger alterable memory 

1 .2 Equipment 

NOVA central processor 

1 .3 External Subroutines 
Unsigned multiply (.MPYU) 

1.4 Other 
None 

2 • Q^^J^a1G_proc?dur^ 

2.1 Calling Sequence 

JSR .DMPY 

address of first word of two word multiplier 

return 

2.2 Input Format 

One double precision operand must be provided in hC0 , 

AC1 (high order, low order) . The address of the 
first word of the second operand must be supplied in 
the word following the JSR .DMPY. 

2.3 Output Format 

The quadruple precision product is stored in four 
consecutive memory locations, highest order to lowest 
order. AC2 will contain the address of the first, 
and highest order, word of the result. 

2.4 Error Returns 
None 
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2 . 5 S tate of Ac tive Registers upon Exi t 

All accumulators and Carry are destroyed* 

2.6 Cautions to User 
Hone 

3 - DISCUSSION 

3,1 Algorithms 

The sign of the result is initially determined 
according to the algebraic rules for multiplication 
of signed numbers. The absolute values of the 
operands are used to determine the absolute value of 
the result. If necessary, this result is negated. 

Assume the absolute values of the operands ares 

2**16 * A + B 

and 2**16 * C + D 

The product of these operands is; 

2**32 * A * C + 2**16 * B'*'C + 2**16 * A *'D + B •* D 

where the single precision products can be formed 

using the unsigned multiply routine. Thus four 
single precision multiplies and three additions are 
required to compute the result. 

3 . 2 r.i rei.tations and Accur acy 

The routine is exact. 

3.3 Size and Timing 

Double precision multiply and double precision divide 
are supplied as one routine. The size for both 
routines is 213 (octal) words. 

The average execution time is 257 JU seconds in addition 
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to four unsigned multiplies at 340 U seconds each. 
Total average time is, therefore, 1.62 milliseconds. 

3-4 References 

Write-up 093-000011 describes the double precision 
divide subroutine which is provided along with .DMPY, 

Write-up 093-000015 describes the unsigned multiply 
routine . 

3.5 Flow Diagrams 

Not applicable 

4. EXAMPLES AND APPLICATIONS 

An ASCII source tape (090-000016) of .DMPY is provided 
with the NOVA software. This tape should be edited into 
a user routine that requires double precision multipli- 
cation. 

5 « PROGRAM LISTING 

A listing of .DMPY follows. No origin has been given, 
enabling the user to edit this routine anywhere into his 

program. 



J 316«£tJ UCJj'rj;^ ."^Li'lP^ 

I "VUZPHtS fW SKi.Vi-O* LfJ',-£LE ^-JCiSUN NUMBERS 

J INPUT: HI IN AC 3 ^JjH *;*>E*. ;* Av I CLJ^ '.MJErt) 

I ii> r.;Ir, i£»j fO dY toOKU AF TEri J5« •UMr'Y 



i output? 
i 



-J : * J 2 l ft 5 f *J k A 62 P 1 N r £ -J I d Y AC 2 
AC2 POINTS TO HIGHEST OrtDEri 



I CALLING SEQUENCE* 

I 33K •i/f.-'Y 

I Al>t)«£33 Of 02 

I HEEUiiM 



I 0E3f«0¥£tJt 



ALL ACS AND GAKrtY 



I KEQJIKESl 



• MPY'J imSIUmO MULTIPLY > 



#8001 

8gg§4 

80006 
00S0? 
008 li 
00811 
#0012 
0§013 
g00 14 

0001s 

0001? 
00020 
00821 
00g22 
fe§0023 
@i#24 



£S>4076 

152400 

04007? 
§44100 
§34076 
010076 
035400 
§2140fe§ 
025401 
004045 
!5!£?,0 

l vHi ;a 

0501 10 

040! a I 
044 - ^2 
030103 
00SU3 
844100 
024101 
0361 14 



00025 030877 
§0026 04»077 
Si§27 121030 



00030 

00i3I 
0«832 
80833 
i§034 

i§0§36 
i@g37 
00040 

00041 
80i#42 
§3043 
80044 



02-6 122 

0$6i;4 

024101 
006114 

830077 
147022 
101420 

i3f#l02 
6134 li30 
'4 ! 4 1 | 

000066 



OMPYf STA 3^.6C03 
JSM *tid50 

STA 1*«£C10+1 

LDA 3*«tSC03 
25d .stCfeI3 
LDA 3#0j»3 
LDA 0*0*3 

LOA I, I s3 
35K *Bd50 

MOVK gj»2 
SU8L 2#2 

3TA 2#*BG13 



Si' A 
STA 
LOA 

3S« 
Si A 

LOA 
3S« 



0#.bC10*2 

1 # «dC10*3 

««c3G3# 
l»«bC10+l 

3 *dC3i 



LlJft 2#»«Cl«d 
STA 0#.BC1# 
MOV 1#0 



LHA 
3SK 
3TA 
LUA 
33K 



i#.tJCI0*3 

0.BC31 

I* .felGlfe}*3 

l#«fc*C10+2 

e*fc<C3i 



LOA 2#.BC10 
Al>t3g 2#1#S2G 
INCiS 0*0 
UXA 2* *bC»kJ*3 
LOA 3#»tfC10^1 
Q$6 .6*013 
JMP »B853 
JMP *BB54 



f save ketukn 

I CLEAK RESULT FLAG WOKO 

I FOKM ABSOLUTE VALUE Of 01 

I SAVE A8SC0I ) 

I KESTOriE AC3 

1 BU HP PAST AOiMESS CONSTANT 

I GET AD0ME35 OF 02 

I GET 01 

1 FOKM ABSOLUTE VALUE Of 02 



1 FOK RESULT MINUS* 1 FQH 

I PLUS 

I SAVE ABSCH2) 

J 01 CLOW> * 02CL0W) 

J UNS16NE0 MULTIPLY 

I SfOHE LOWEST QHOEd OF KESULT 

I Ok 'lUV) * i)2<HiGH> 

I ADD HIGH QrtDEK TO FIRST CM035 

I PttODUCT 

; OiCHIGH) 

I 3AVE HiDH ORDEK OF 1ST CHOSS 

J LOW OHUEit TO BE AOOED TO 

I 3EC0M0 CM033 

J U2CLU^) * 31 CHI3H) 

I SECOND CK033 PKODUCT 

I 2N0 LOWEST OHOEK OF KESULT 

I 0KH1GH) * D2<riiGH) 

I AOO M16H OKtlllM OF CK0S3 TO 

J rflGH OiilJEM MULTIPLY 

I AOO HIGH OiiOEH OF 1ST CHOSS 



I GET LOto QHiJEH RESULTS 

I TEST 31G^ OF KESULT 

I NE3ATIVS. - NEGATE rtESULT 

1 POSITIVE - 5TCWE RESULTS 



I £NTf<y hi *r^-0 FOrtMS ABSOLUTE VALUE Of NUM8EK 

I IM ACt/. AC! A-mj; INCREMENTS AC 2 

I ENTrfY AT #8856 NEGATES AC0# AC1 AMD INCREMENTS AC2 

I ENTKY AT .BBS! NEGATES AC§# AC I 



£C5??*i I, %S t S 1 3 ♦•:?£;? £ ? MOVL^ 0#0*SNC 

0S?-45 f:*:"::^//^ jr*P @*3 

S#34 - U : * Z ' * =itf5 f ! I r: C * p 2 

t ? ?^5 ::• ' ?A t?4 * 3r-i is N£ G 1 # 1 * iitSK 

00051 100001 COM 0#0*SKP 

§§§S2 13§4§§ MEG 2*0 



I ENTKY AT .8BS2 FOKMS ABSOLUTE VALUE OF NUM9EH 
I IN AC#* AC1* AC2^ AC3 

I ENT«Y AT .BBS3 FOK^S 0UA0KUPLE NEGATE 

1 ENTMY AT .BBS>4 STOKES THE ACS IN A BLOCK STAKTIN6 
I Af *8Ct0 

I IF *DMPY CALLED TME ROUTINE* CAMKY WILL BE 2EHO 

i 

I IF .001 V CALLEO THE KOUTINE* CAt^HY WILL BE N0N-2EH0 



#«§S4 


101.113 


*88S2 8 


MO¥L# IS#0#SNC 


§§§§§ 


i§§§66 




JMP .8854 


#3§56 


1 7^4-34 


•B8S3S 


NEG 3r3*S*SK 


§§§5? 


1£??0I 




QOr. 2#2^SKP 


00060 


15^464 




NS6C 2*2j»S2rt 


£536 1 


12*' •?■''?. 




ca;*i i*i*s-<P' 


0#§62 


1244 64 




NEGC 1»1*S£K 


§§§63 


i^e-?si 




COf* i#i*SKP 


§ncf.4 


i3?*6S 




NEGC 0#i#SNM 


g0§65 


101060 




MO¥G @#i 


0f§66 


e4g§7? 


*BB54f 


si a 0#.etn.? 


§9067 


0441 Si 




ST A i##BC13+l 


'•"'"■f.?^ 


®501il 




STA 2*»SC10+2 


c:r7?. 


054 102 




STA 3»*BC10*3 


i§e?2 


101282 




MOV §j»ij»Si-C 


0#«?3 


000136 




J«P «BC98 


§8074 


030115 


•BB55I 


LyA 2#.8C32 



0§075 §§2§?6 



JMP i.BC03 



I rtETUKN TO DIVIDE KOUTINE 

J POINTEK TO MULTIPLY RESULTS 

I QH DIVIDE KEMAINOEft 

I ME TURN 



000 76 000000 .BC038 

000004 *«C10! *dLK 4 



000002 .dCllS •BLH 2 
0:j 0.300 .eC!&5 .8LK 3 
00110 000000 .BCI3S 



Hi 003000 • B C 1 4 3 

00112 00004! .8020 S 41 

00113 000214 -BC30S »MPYU 

00114 000215, .BC31S ««PYA 

00115 000077* .BC32I .BC10 



I SAVE MEfUHN 

J SAVE Ado<U>* 

I OM KEMAINDEK 0'*' OZVICE. 

I SAVE AriSC V> 

J fEMPOKAKY STORAGE 

J i 1, 'itf S \X- J U - / Ki£ ^A 1 Ni>£ K 

J FLAGS* 

J OH SIGN Of *wtfIPLY rtESULF 

J 1,[£K A r ! ^4 C Oil N r WOK 

I IIEKAflON FOK DIVIDE <33> 



J UNSIGNED MULT 1 PLY 

1 UNSIGNED MULTIPLY AMD ADO 

J PQiNFEK ft) KESULTS 



♦ft 



SIGNED uQ'JX'-Z £ IV JOE 

£* i v : c:w : : - a s i ? * »£.; , */jlt : ^12 - e hcc i s i on numsers 

INPUT* DIVI-Cc^ <<i>* A QUADRUPLE PRECZ3IQ.N 

POINTS ro ihZCit* QH^ii ¥QR$} BY AC2 

A". J OCCUPYING FOUR CONTIGUOUS MEMORY LOCATIONS 

DIV'SCI* {V>> A. DOUBLE PESO'S ION TWO f S 
COMPLEMENT NUMBER IN ACi# AC 1 <KIGW# LOW) 

OUTPUT s :M pQ 1 1 U2#U3/V3* VI 

DOUBLE PRECIS* 0-1.- SXSMS3 QUOTIENT IN AC0# AC 1 
CKI-3H* LOW) 

DOUBLE PRECISION 516WE0 CSAME AS DIVIDEND) 
REMAINDER POINTED* TO BY AC2 <P;.GK# LOW) 

CALLING SEQUENCES 

JSM .DO IV 

RETURN 



EHHOti CONDITIONS! 



•If A8S€U«3#U1> >= ABS<V0*V1)« 

CARRY WILL BE 

SET AND RESULTS UNPREDICTABLE 

IF THE QUOTIENT MAS AN ABSOLUTE 
VALUE >« 2**31 » 
CARRY WILL BE SET AND 
RESULTS UNPREDICTABLE- 



J DESTROYED* 



ALL ACS AND CARRY 



m> 16 054016 .DDIVS 
001 17 0SS077 
§§128 IS8S20 
0:3:21 034045 

00122 #4il»3 .BC99t 
#0183 044104 
g§124 I5SI20 



S»I25 


03SS77 


05.* 26 


i2l0i« 


#i!27 


101112 


00130 


1 75400 


031 32 


0541 10 


00132 


S25i§l 


00 1 33 




8§134 


@31g§2 


B§t3S 


§§i@S4 



STA 3#»BCe3 


1 


STA 2#*BCii 


1 


su?-;i^ 2#a 


* 


JSR *BB5i 


1 




. 1 


STA 0»«BCU 


1 


STA 1#»BC11*1 




MOV2L 2#3 


1 



LDA 


2#.scia 


LOA 


2*0*2 


m¥U* :-0*S-5C 


XNC 


3#3 


STA 


3#.BC13 


LPA 


1*1*2 


LOA 


3*3*2 


LDA 


2*2*2 


JMP 


• 8852 



SAVE RETURN 

SAVE AC2 TEMPORARILY 

CLEAR FOR RESULT SIGN FLASS 

FORM ABSOLUTE VALUE OF 

DIVISOR 

smz ABSCV) 



SIGN TO BIT 



POSITION DIVISOR 

# 14* SET CARRY 
I RESTORE AC 2 

I GST DIVIDEND 



I SIGN FLASS ARE SET 



I FORM ABSOLUTE VALUE OF 

I DIVIDEND (NOTE CARRY SET > 



., 4 . , 







id id I 

00 i 

00 i 

001 

001 

001 

id Id I 



103 

000145 
142.333 

166432 
002076 
•:*•'..•* '» ; 2 
«S 4r--i 1 1 

4? 101020 
30 000 1 66 



4k3 



3-</-^ '■ >i * '■•■■ « c* '^3 

£ J r * I - 2 .» 2 / 3 xu 

JflP 8.BC03 
LOA 3««dC20 

3 < / -- ii •* - = 3 i '•'• 
MOV 2 0*0 

3MP .BG90 



W6\o\ 166423 .BC89J SUB if 3#t#5NC 



00152 142021 

00153 142420 

00154 008166 

00155 125100 

00156 101100 

00157 030103 

00160 034104 

00161 142413 



00162 
00183 
00164 
00165 
00166 
00167 
00170 
00171 
00172 
00173 
00174 
00175 

00176 
001 77 

00200 

00H 01 
00202 
00203 
00204 
00205 
00206 

00207 

00210 
0021 1 
00212 



000166 

i^2^?i3 

166412 
000151 
034102 
030101 
175100 
151 100 
050101 
054102 
0141 1 i 
000153 

0301 10 
151202 
00404 7 
040077 

0«'«:00 

020101 

1 13102 
002076 
131202 
004050 
I0I020 
0000 74 



5c36 2^0 

J«P .8C90 

• acms movl i#i 

MOVL 0*0 

LDA 2#*BGi 1 
LOA 3,**3CI 1*1 

sua# £>0»snc 

JMP .8090 
SUfc*# 2*0;»SNK 
SUtM 3#i#S2C 
JMP *BC89 
■ BG901 LDA 3#«*k§Gie + 3 
LDA 2#»BC10+2 
MOVL 3 #3 
MOVL 2*2 
sTA 2-- *rf>10+2 
SFA 3#*BCI0*3 
03* -BC14 
JMr» *BC88 

LOA 2#.BC13 
MOVK 2*2,r,^C 
J3K .fc«5fe 
i«TA 0# »b0 1 
a f A I * •SG10+1 
LOA 0# s fcJCl 0*2 
L0A l*»bCi0+3 
MOVL 0#3#SiC 
JrtF t»BC03 
MOVK 2* 2-«3 A 3 

^J7.< 0*0 



I %e^ j;v 3 'jr.- -- jk EKKQH TEST 



1 EKrtOK* CArirtf IS SET 

I +33 

I STOKE LOOP COUNT 

I CtEftK CAHhY 



I SUBFiiACr Oi^IrW, FKO^ HIGH 
I DIVIDEND 



I GET 01VI3OK 

I WILL DIVISOR 60 INTO 

I D1VIOEN0 

I NO 



I YES 



1 SHIFT DIVIDEND AND QUOTIENT 



IMP 



• BfciSS 



1 NOT 0ONE VET 

I GET FLAGS 

I NEGATE *EMAIND£rt# WC AC2 

I AND STOrtE If 

1 GET LVOViE^-:" 

I TEST FOrt >2**3i 

1 YES* CArtKY SET FOri £HriQK 

I NEGATE (iUOTIE^r 

I CLEAN GAKHY 

1 KETUKN 



